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Abstract 

We present results for the fermion determinant in an external Yukawa-coupled 
scalar field, a quantity relevant to instanton and sphaleron transition rates as well 
as bubble nucleation occuring in the electroweak phase transition when the scalar 
field is identified with the Higgs field. We calculate the determinant exactly and 
find an approximation (essentialy of the gradient expansion type) that reproduces 
the exact results very well if the fermion mass irtp becomes larger than the inverse 
of the size R of the background field configuration. The results are given for 2, 3 
and 4 Euclidean dimensions. 
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1 Introduction 



The fermion determinant in an external Higgs field is a quantity that enters 
semiclassical rates [|I|, [| for sphaleron |3], |], [5] and instanton || transitions as well 
as for the nucleation of bubbles or droplets f?|, |8| in first order phase transitions, 
as possibly the electroweak phase transition occuring in the early universe. The 
classical action is usually minimized by spherically symmetric field configurations, 
the mass term produced by the Higgs field is then given by a radial profile 4>(x) = 
f(r). We will consider here field configurations that approach the broken Higgs 
phase as r — > oo. 

The exact calculation of such determinants has to be performed numerically 
since usually even the background field configuration is only given as a numerical 
profile. A very efficient method based on a theorem on functional determinants 
has been developed recently by V. G. Kiselev and one of the authors (J.B.) 
fLOfl . We will apply this method here for evaluating the fermion determinant 



exactly in 2, 3 and 4 Euclidean dimensions. 

It is certainly useful to dispose of analytical approximations to the fermion 
determinant, especially if the back reaction of this determinant on the dynamics 
of the Higgs field is being considered in the form of an effective action. It is 
our main purpose here to present such an approximation and, using our exact 
results, to demonstrate that it approximates the fermion determinant very well 
for sufficiently large fermion masses. This covers indeed the most relevant case 
since the effective action induced by fermion background fields behaves essentially 
as nip where D is the Euclidean dimension and so is most relevant for heavy 
fermions. 

The paper is organized as follows: In section 2 we recall some general rela- 
tions and definitions. In section 3 we present the spherical decomposition of the 
relevant Dirac equation in D = 2, 3 and 4 and discuss the numerical evaluation 
of the determinant. In section 4 we develop our approximation. The numerical 
and approximate analytical results are presented and discussed in section 5. De- 
tails on the Dirac operators and the partial wave decomposition are given in the 
Appendix. 



2 Basic Relations 

The Dirac operator for a fermion with Yukawa coupling gy to an external scalar 
field takes, in Euclidean space, the form 

V = 7^ - m F <j)(x) = P (0) - m F {(j>{x) - 1) (2.1) 

where we have normalized the "Higgs" field to its vacuum expectation value v 
so that rrtp = gyv and where T>^ is the free Dirac operator in the "broken 
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Higgs phase" where (f>(x) = 1. The gamma matrices are hermitean and satisfy 
7 M 7„ + 7^7^ = 25 fiu . We also define the positive definite operator 

M = V^V = -d 2 + m 2 F + V(x) (2.2) 

where the potential or vertex matrix V takes the form 

V(x) = m|(0(x) 2 - 1) + m F ^d^(x) (2.3) 

We can evaluate the determinant of the Dirac operator, normalized to the free 
Dirac operator in the form 

lndet(P/P(°>) = | In J (2.4) 
-d 2 +m 2 F + V(x)\ 



J = det 



-d 2 + 



m 
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Furthermore the fermionic fluctuation determinant is related to the corresponding 
one-loop effective action by the relation 

S eff [4>] = - l -\nJ (2.5) 

The case D = 2 could be relevant for lattice simulations of bubble nucleation 
at low temperature. The case D = 3 will be considered with hindsight of an 
application as high temperature limit of the four dimensional theory. The four 
dimensional case which would be relevant to bubble formation at low temper- 
atures (instanton type transition). The details on the 7 matrices are given in 
Appendix A. 



3 Numerical Computation of the Fluctuation 
Determinant 

For spherically symmetric background fields <f)(x) = f(r) the determinant decom- 
poses into contributions of partial waves in the form 

In J = Y / dx(D)\nJ x (3.1) 

where A denotes a coupled channel of partial waves, J\ is the partial wave deter- 
minant and d\(D) is the degeneracy which depends on the dimension D. J\ is 
related to the partial wave reduction of the operator M. which takes the general 
form 

M A = Mf +V(r) 
» .(o) d 2 D - 1 d L? 2 

= -53-— * + F +m ' (3 ' 2) 
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via 

JX = ^ (m|) (3 - 3) 

The degeneracies d\, and the matrices specifying the centrifugal barrier and 
V(r) are given in Appendix A. 

For the numerical computation of the determinant J\ we use the method 
described in Ref ||10|| , generalized to a coupled channel system. It is based on 



a theorem on functional determinants 01 which we recall below. The theorem 



has been used for a similar purpose for the case of massless fermions . In the 
following we consider the contribution of one channel A. 

Let ipi a (r) be a fundamental system solutions of the differential equation 

M x ip ia = (3.4) 

and let ip^ be their free counterparts, both regular at r = 0. Since there are 
as many independent solutions labelled by a as there are components labelled 
by i these fundamental systems form quadratic matrices denoted by \& = {ipi a } 
and \E ,< -°' ) = {^-^ } respectively. Let furthermore the solutions be normalized such 
that 

Um*(r)(* (0) (r)) -1 = 1. (3.5) 

Then according to the theorem on functional determinants || the partial wave 
determinant is obtained as 

d6t (Sif ) = ™ det^°)(r) (3 - 6) 

In the cases considered here — the Dirac operator with a scalar mass term 
in two, three and four dimensions — the partial wave reduction results in simple 
(2x2) systems. The "free" differential operator M A °^ is diagonal and the solutions 
are modified Bessel functions which be denote as b\ for the upper component; 
their detailed form is given in Appendix A (see Eq. ( A..1C)). The solution for the 



lower component ist &a+i in all cases considered here. 

The solutions of the full operator Ma are now written as 

i>ia{r) = bx + i-i{r){8icL + h ia {r)) (3.7) 

where the boundary condition at r = requires the functions hi a (r) to vanish 
there. Since the Bessel functions can be taken out of the rows of the matrix ip 
when forming the determinant they cancel in J\ and we have simply 

J A = lim det(l + h(r)). (3.8) 
As in [irD| we write directly a differential equation for the functions h ia . It 



takes the form 

A A h = V A (l + h) (3.9) 
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where the differential operator Aa and the modified potential Va are defined in 
Eqs. ( |A. 18| - |AT20|) . This differential equation is inhomogeneous and h may be 



expanded w.r.t. orders in the potential V. Denoting the contribution of order k 
in the potential as h^ fc ^ we introduce also the contributions summed from k to oo 

as 

oo 

h (fc) = £ h (0 ( 310 ) 



l=k 



so that h = h^. The differential equation can then be used to get the functions 
h of various orders recursively via 

A„h (fe+1) = V„h (fc) 



A n h (fe+1) = V„h (fc) (3.11) 

where further h(°) = 1. Using these contributions of specific order in V it is easy 
to get rid of the leading, divergent contributions to In J\ and therefore to In J . 



This procedure has been discussed in [10] and can be performed analogously for 



a coupled system: In D = 2 and D = 3 the contribution of order k — 1 has to be 
omitted, in D = 4 the orders k = 1 and 2; the finite parts of lnJ\ are then given 

h y _ 

It J A ) (2) = lim (Tr[ln(l + h) - h] + Tr(h (2) )) (3.12) 

and 

(In J x )^ = lim ( Tr[ln(l + h) - h + -h 2 ] 

1 — >oo \ 2 

+ Tr(h^) - ^Tr(h^(h« +h«))) (3.13) 

The matrices in square brackets can be diagonalized in terms of the eigenval- 
ues of h and reduce then to "subtracted logarithms" which can be evaluated 
with arbitrary precision. The remaining terms are separately of order (2) and of 
order(3) in the first and second equation respectively. These expressions do not 
contain therefore any dangereous small differences of large numbers. The total 
fluctuation determinants are denoted correspondingly as {hij)^ and (InJ')^. 

The divergent contributions have to be regularized and renormalized and this 
depends on the specific model under consideration. These contributions can 
easily be calculated analytically, at most one needs the Fourier transform of the 
profile (see e.g. |1C| where this step was included). Here we are interested in the 



comparison with the large mp expansion and for this purpose it is sufficient to 
consider those finite parts of the fluctuation determinant which are obtained by 
omitting the renormalization parts completely. 

The numerical computation has been carried out with all the technical details 
as described in Ref. ]ll|. The results will be discussed in section 5. 
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4 The Large Mass Approximation 



If the fermion mass mp becomes much larger than the inverse size R~ x of the 
background field configuration one expects the Dirac operator to "seize" the 
background field only locally. Therefore it should be possible to expand the 
fluctuation determinant with respect to derivatives which on dimensional grounds 
results in an expansion with respect to 1/mp. Various techniques are available. 
While the heat kernel expansion is used frequently we use the ordinary Feynman 
graph or resolvent expansion here. While heat kernel expansions are often broken 
off at a fixed power of the heat kernel time we will keep track here of all terms 
of the leading and next-to- leading power of the mass mp. 
The fermion determinant can be written exactly as 

ln l 7 = Trln(l + G (0) V) (4.1) 
where the free Green function is defined as 

MG {0) = 1 (4.2) 
We define the Fourier transform of the potential V as 

V(q) = J d D xV(x)exp(-iqx). (4.3) 
The fluctuation determinant can then be expanded as 

where 

Q; = X> (4-5) 
i=i 

The range of the integration in the variables is determined by the inverse range 
of the external field configuration R~ l . If mp ^> R^ 1 the scale for the integration 
in p is set by the fermion mass mp. This integration yields therefore a factor 
mp . Recalling the form of the potential 

V = m^(0 2 - 1) + m F7 <90 (4.6) 

we see that in leading order the product over the potentials yields a factor m 2 p 
so that the leading behaviour of YnJ is m® Q The next-to-leading order term 
is obtained by taking into account the first nonleading terms in the expansion 

1 Note that this power counting is different in the case of chiral mass terms of the form 
mo exp(i75</> a r a ) where only the gradient term linear in mo appears in V. 
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of the denominators and of the trace of the product over the potentials. Using 
symmetric integration we have 



n 1 1 r n Q 2 

+ ±v£&(L<** + ±*j\ (4.7) 

up to terms of order m F 2n ^ . Similarly we have 



TrlpM ~ ^[n(0 2 -l)(q,O 



J 3 

- XYs^i'Mmm n (^ 2 ~!)M (4.8) 

F l>l> 

up to terms of order m^ n_4 ' ) . Here vp is 2 for D = 2 and 4 in D = 3 and D = 4. 
Inserting these expansions one obtains after transforming back to x space 

In J = vpm^Y,- " / d D xUl{D,0,n){<f) 2 - l) n 

n=l U ^ 



+ ^^/p,O,n)(0 2 -ir 2 (^0) 2 (4.9) 



"F 



2{n-\)n{n - 1) , ; , n , , 2( ,v<-2, ;) ^2 
3m 2 



-J(AO,n + l)0 2 (0 2 -ir 2 (^) 5 



+ 



where 



I(D,m,n) = / 7 xn/i r ; . (4.10) 
V ' ' ; J (2n) D (l 2 + 1)« V ; 

The sum over n can be performed after evaluating explicitly the (dimen- 
sionless) integrals I(D,m,n) in the various dimensions. The first few terms are 
however formal only: the integrals are divergent. As for the exact computation 
we have to omit, for D = 2 and D = 3 the order n — 1 and for D = 4 the orders 
n — 1 and 2. We denote the summation over n = n , ..oo by the symbol (n ) as 
already in the previous section. 

Writing the result in the form of an effective action we have for D = 2 



2 



Sf^] = - l -\nJ^ c ^/^(l + 2 (ln(0 2 )-l)) 

24tt 



^ (4-11) 
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for D = 3 



^3 X (M! (412) 
127ri |(/)| V ; 



and for D = 4 



^ 3 / ) /M=4 ln ^ (3) " -^/^(0 4 ln(0 2 )-|(0 2 -l) 2 -(0 2 -l)) 

|rf 4 x(^0) 2 (ln(0 2 )-^(0 2 -l)) (4.13) 
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5 Results and Conclusions 

We have evaluated the exact fluctuation determinants described in the previous 
sections in D = 2, 3 and 4 for two different typical profiles, a "small" and a "large" 
bubble. For convenience we used an analytical expression for the radial profile 

f(r) = (1 + exp(-Vr 2 + 1 + VR 2 + 1)) _1 (5.1) 

This function behaves near r = as C\ + C2?" 2 in analogy to the solutions of the 
typical field equations for the Higgs field, for r — > oo it behaves as 1 — C3 exp(— r) 
as it would behave for a Higgs field of mass 1. R determines the point f(r) = 1/2 
and therefore the typical size of the profile. We have chosen R = 2 and R = 6. In 
the second case the function f(r) goes to .006 at r = so that the symmetric Higgs 
phase is almost reached. This is a situation where breaking off the expansion in 
(0 2 — 1) at a finite power of this quantity would lead to a poor approximation. 
The profiles are displayed in Fig. 1. 

The results of the exact and approximate computations are presented in Figs. 
2-7, both for the small and large sized profiles and for D=2, 3 and 4. 

The precision of the exact computation can be estimated from the convergence 
of the partial wave expansion, from varying the grid density for the radial variable 
(we have typically 1000 points in the range < r < 15) etc. to be better 
than of 10~ 3 . The agreement with the estimate at large fermion masses further 
corroborates this estimate. The precision becomes poorer UttifR exceeds a value 
of ~ 12 since then the asymptotic regime of the partial wave series, requiring 
I ^> mpR is not yet reached at our maximal value of / = 25; while the higher 
partial waves are included using a fit a// 3 + b/l A the parameters of this fit are 
then not reliable. 

In the case D = 2 and R = 6 a bound state eigenvalue approaches in the 
S-wave and the calculation becomes numerically delicate at mp > 3. The poor 
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agreement in next-to-leading order between asymptotic estimate and exact results 
which is apparent in Fig. 3 is however not due to the numerics which we have 
checked in all aspects. The <fr~ 2 singularity in the kinetic term of the effective 
action for D = 2, Eq. ( 4.11 ), indicates a bad convergence of the 1/mp expansion 
if cf) becomes small, even if, as on the present case, the integral exists. In all 
other cases, for both profiles and all dimensions, the first two terms of the l/mp 
expansion agree very well with the exact results for sufficiently large vrtp. The 
asymptotic regime sets in earliest for D = 4. 

In conclusion we have found in general an excellent agreement between our 
exact and approximate evaluations of the fluctuation determinant of a heavy 
fermion coupled to an external Higgs field, the singularities of the gradient terms 
at 4>(x) = are however potentially dangereous and need further consideration. 
The comparison corroborates the precision of the numerical method. The analytic 
approximation has the form of an effective action including all terms with at most 
two derivatives in the external field. Its should be interesting to include it into 
the action from which the classical profiles are determined. 
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A Appendix 



We shall summarize the notations and formulas in D=2,3 and 4. The Euclidean 
Dirac operator, V, is defined as 

T> = 7 A - m F (j)(x) (A.l) 
As gamma matrix convention for D = 2 we take 

71 = °~i 

72 = 0- 2 

75 = -o"3 

llJ.lv = - ^^75. 



For D = 3 and D = 4 we use 



7, = ( s °* t ) (A.2) 



2 The three dimensional theory is considered here as a three dimensional reduction of a four 
dimensional one as it occurs in finite temperature quantum field theory 
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with 



E„ = j 12x2 ^ " 4 1 „ (A.3) 



75 is then 

This is just the chiral representation. For a field configuration which de- 
pends only on the radius, <j){x) = f(r) the "squared" Dirac operator M. defined 
in Eq. (|2.2|) is given by 

M = -d 2 + m 2 F + m 2 F (f(r) 2 - 1) +m F f'(r)x^ lx 

d 2 D-ld L 2 2 w , 
+ — + m| + V(r) 



where x M is the unit vector of the coordinate x M , potential V becomes 

V = m|(/ 2 -l)+m F /%7 M (A.5) 

and L 2 /r 2 is the centrifugal potential. 

We take as a basis for the solutions of the Dirac equation in D = 2 the spinors 

where the angular momentum n takes the values n = 0, ±1, ±2... For D = 3 we 
use the spinors 

*=(1, L ) = (-?V X V (A. 7) 

For the four dimensional theory we take 

-(£)-(*&) (A - 8) 

X is a 2-component spinor multiplied by a spherical harmonic function and 

L 2 X = 4j(j + l)x (A.9) 
where j(= 0, 1/2, 1...) is the magnitude of 'angular momentum' in D=4 and 

L 2 x=](] + l)x (A.10) 

where j(= 0,1,..) is the magnitude of angular momentum in D=3. Then the 
'squared' Dirac equations are 

1 — 23 A 2 — 3 2 
- u" ^—u' + — -f- u + m 2 F u + m 2 F (f 2 - l)u + m F fv = (A.ll) 
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_ v » _ LJ^ V > + (A + 1)2 + m 2 ^ + m 2 (/ 2 _ 1)v + fflpf'u = (A.12) 



where 



' \n\ = 0,1,2,.., for D = 2 

A = i j + 1/2 = 1/2,3/2,5/2,.., for D = 3 (A.13) 
k 2j=0,l,2,.., forD = 4 



2 - D 

and (3 is — - — . In this basis the potential is simply the (2 x 2) matrix 



V 



m 2 F (f 2 — 1) m F f 
m F f m F (f 2 ~ 1) 



(A. 14) 



For the free equation(/ = 1), the regular solutions are 

*i = 

where 



b x (m F r) \ 

) 


b x+1 {m F r) 



(A.15) 



b x (m F r) = r^I x (m F r) 



(A.16) 



and I x (m F r) is the modified Bessel function of order A. Therefore, the solutions 
of the 'squared' Dirac equation are expressed as 



*Pia = (Sia + h ia (r))b n - 1+i (m F r) 



(A.17) 



where a labels the different solutions and % the components of the solution. The 
functions h ia (r), are determined by 



where 



and 



— + (— + 2m r b '^ {mFr) ) - 
dr 2 \ r F b X -i +i (m F r) J dr 



V\,ij - v if 



b X - 1+j (m F r) 
b\-i+i{m F r) 
The degeneracies of the solutions are given by 

2, 



(A.18) 
(A.19) 
(A.20) 



for D = 2 

d x (D) = { 2(2A + 1), for D = 3 

2(A + l)(A + 2), forD = 4 



(A.21) 
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Figure Captions 



Fig. 1 The Higgs field profiles. We plot the radial profiles given in Eq. ( |5.1|) for 
R = 2 and R = 6 as used in our computations. 

Fig. 2 The fermionic effective action S^/mp for D = 2, R = 2. The solid line 
represents the analytic approximation of Eq. fl4.ll ), the numerical exact 



results are plotted as points (•) connected by straight dotted lines. The 
straight dashed line indicates the limit mp — > oo. 

Fig. 3 The same as Fig. 2 for R = 6. 

Fig. 4 The fermionic effective action S^h/m Z p for D = 3, R = 2. Description as 
for Fig. 2, the analytic approximation being given in Eq. (|4.12j ). 

Fig. 5 The same as Fig. 4 for R = 6. 

Fig. 6 The fermionic effective action S^h/m F for D = 4, R = 2. Description as 
for Fig. 2, the analytic approximation being given in Eq. ( |4.13j ). 

Fig. 7 The same as Fig. 6 for R = 6. 



12 




Figure 1 




Figure 2 



10 




Figure 4 
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